Theory of the spin Hall effect, and its inverse, in a ferromagnetic metal near the Curie 

temperature 
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We give a theory of the inverse spin Hall effect (ISHE) in ferromagnetic metals based on skew 
scattering via collective spin fluctuations. This extends Kondo's theory of the anomalous Hall 
effect (AHE) to include short-range spin-spin correlations. We find a relation between the ISHE 
and the four-spin correlations near the Curie temperature Tc- Such four-spin correlations do not 
contribute to the AHE, which relates to the three-spin correlations. Thus our theory shows an 
essential difference between the AHE and ISHE, providing an essential complement to Kondo's 
classic theory of the AHE in metals. We note the relation to skew-scattering mechanisms based on 
impurity scattering. Our theory can be compared to recent experimental results by Wei et al. [Nat. 
Commun. 3, f058 (2012)] for the ISHE in ferromagnetic alloys. 
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The spin Hall effect (SHE), which converts charge 
current into spin current, along with the inverse effect 
(ISHE) which re-converts the currents, is one of the key 
phenomena for the further development of spintronic de- 
vices [H, . A difficulty in exploiting the effect is that it 
depends on a spin-orbit interaction which, as a relativis- 
tic effect, is intrinsically weak except in heavy elements. 
Recently the doping of gold with low densities of iron 
impurities was proposed as a mechanism to enhance the 
SHE through spin fluctuations of the individual iron mo- 
ments [1, This raises the question as to whether the 
effect might also be sensitive to collective magnetic fluc- 
tuations, in particular, close to a magnetic phase transi- 
tion, where the fluctuations are strong, in analogy with 
the anomalous Hall effect (AHE) in a pure ferromagnetic 
metal such as Fe Q and Ni Q , where a peak of the Hall 
resistivity appears below the Curie temperature Tc- Wei 
et al. Q observed a feature in the ISHE, scaling with T c 
for different samples of a metallic alloy, which suggested 
that this is in fact so, but the form is quite different from 
that of the AHE. 

The observations of the AHE near the Curie temper- 
ature were explained by Kondo in his calculations [8| of 
an s-d (or s-f) lattice where, apart from the standard 
Coulomb interactions between the orbitals, he included 
on-site spin-orbit interactions. This gives rise to skew 
scattering in the conducting s band. We will show here 
that while both the AHE and SHE are generated by skew 
scattering that is linear in the spin-orbit term, they have 
qualitatively different behaviors near Tc as they depend 
on correlations of distinct order. In Kondo's formula- 
tion Coulomb interactions are included by keeping s- 
and p-scattering channels for the conducting electrons 
around each localized orbital. The angular momentum 
can thereby be transferred from the local moment to the 
conduction electrons. As we are interested in the effects 
of skew scattering, we will write only terms depending 
on the vector product re' x k, where re (respectively re') 
is the unit vector in the direction of k (respectively k'), 



the momentum of an electron in the s band. The s-d ex- 
change interaction is then 0] (omitting terms depending 
on quadrupole moments that do not contribute here) 



H 



(3i/4)F 2 L n • (re' x re) 



+2{S n ■ s c ){ J(k,k') + (3z/2)c 2 F 2 L„ • (re' x re)} 



L n and S n are the total orbital and spin angular mo- 
menta of the localized d orbital on atomic site n at posi- 
tion R„, of a total number of N magnetic atoms in the 
crystal. s c is the spin operator for the conduction elec- 
trons. J(k, k') = Fq + 2T\(re ■ re') with ^0,1,2 exchange 
terms generated by Coulomb interactions between the lo- 
calized d electrons and the different channels of scattering 
with propagating s orbitals. The coefficients c 2 appear 
in the summation of individual orbital and spin angular 
momenta of each electron within the localized orbitals to 
give the total momentum Q. Assuming an orbital sin- 
glet but adding a small spin-orbit Hls = ^^2, n L n • S n , 
it generates operators L n on each site and gives a cor- 
rected Hamiltonian with the coefficient Ai linear in the 
spin-orbit coupling A and determined by the local crystal 
field levels M, 



H = - ]T AT-V^-^a^v 



2(5 n -s c )J(k,k') 



n,k,k' 

+ iA 1 F 2 {S n ■ (re' x re) + 2c 2 (S n ■ s c ) (S n ■ (re' x re)) 
+ 2c 2 (5„ • (re' x re)) (S n ■ s c ) - 4c 2 /3 (S n ■ S n ) (s c ■ (re' x re))} 

(1) 

Considering only elastic scattering, the matrix element 



2 



of the above Hamiltonian is 



k±,k'± 



±J(k,k')(M„-(M n )) 



N-l 



2c 2 



,i(k'-k)-(R„-R ; 



((M n -(M n )) 2 (M 2 ~(M 2 ))), 



+ iAiFaC/c' x /c) c {±2c 2 (Af 2 - (M 2 )) + (M„ - (M n ))} 

(2) 

M ra is the magnetization of the localized electron of the 
n-th ion, and £ is the magnetization direction below Tq 
or the direction of polarization of the spin current above, 
a direction which will be determined by the spin injector 
in the devices used. Subtractions of the lattice-averaged 
powers of the magnetization at each site, (M n ) and (M%), 
appear because any modification of the crystalline poten- 
tial simply rcnormalizes the Bloch functions of the con- 
duction electrons and does not contribute to scattering. 
There are two parts in Eq. ([2]): One is the J(k, k') term, 
coming from the s-d scattering, and the other is the Ai 
term, which is linear in spin-orbit coupling A. ± are spin 
states of the conduction electron, (k' x k)^ comes from 
the approximation of including elastic scattering only. 

In order to find a net current transverse to the volt- 
age, the transition probabilities must be calculated to at 
least the second Born approximation , i.e., the transi- 
tion probabilities are to third power in the Hamiltonian. 
There is a triple summation over sites which are reduced, 
with Kondo's assumption of uncorrelated fluctuations, to 
a single summation. Now we generalize to include cor- 
relations but for simplicity we include those with only 
two-site indices, but all powers of spin. The transition 
probability from k'± to k± is given by 

W{k'±, k±) = 6(E k - E k , ){Ux + U 3 + V ± (k', k)}, 

where the s-d scattering contribution of U\ is given by 

Ui = (2ir/hN)\J(k,k')\ 2 Pl , 



N-l 



(k'-k)-(R„-R ; 



<(M n - (M„))(M - (M ))}. 



N-l 
n=0 

A contribution U3 = (2ir /fiJV)P(k, k') is included to rep- 
resent non-magnetic contributions, static impurity scat- 
tering, or phonons, where P(k,k') is assumed to depend 
only on the angle between k and k , and not to show 
significant temperature dependence close to the mag- 
netic phase transition. The spin-orbit contribution of 
V A± (k',k) is obtained as 

V^.k) = (2ir/h)(A 1 F 2 V/3TrN 2 )(2m 3 ) 1 / 2 h- 3 E 1 k /2 
x (k' x K ) c [(r x ± r 2Q ){3F 2 + 4F 2 (k •«')} + (r x ± r 2b ) 



{-4FbFi -8F 2 («•«')} 



(3) 



where r 1 , r 2a , and r 2b depend on k' — k and temperature, 



N-l 



ri = e i; ( k '- k HR"-^)((Af„ - (M n )) 2 (M - (Mo))), 



'2b 



2c 2 J2 



= i(k'-k)-(R„-R ; 



n=0 



((M n - (M n ))(M 2 - (M 2 n )) 
x(M Q -(M ))). 



These generalize the on-site n = terms, which we denote 
as r° or r 2 , of Kondo. There are now two distinct contri- 
butions at fourth order in spin, r 2a and r 2b . All off-site 
coefficients (n 7^ 0) depend on the momentum transfer, 
making integrals slightly more complex than Kondo's. 

Inverse spin Hall effect. Here there are two spin po- 
larizations possible: first, that of the spin current, in- 
troduced by a spin injector, assumed to have a much 
higher Curie temperature than the Tq of the ferromag- 
netic metal (spin detector) in which the skew scattering 
and the inverse spin Hall effect occur. It is the injector 
which defines the direction of quantization of the spin 
current £ both above and below Tq. The second is the 
ordered moment of the ferromagnetic metal (spin detec- 
tor), which we assume to be parallel to the spin injector 
below To The incident spin currents are driven by the 
spin diffusion field: the difference in gradients of two elec- 
trochemical potentials Ep d, [Io|, T = — -V(ef ± 5ef) 
with magnitude T and direction ±<x for the spin states ±. 
The spin current is scattered to create a transverse volt- 
age. Solving the Boltzmann equations for distributions 
over the two Fermi surfaces of spin-up and spin-down 
electrons using the skew scattering probabilities Eq. ([3]) , 
the current densities H can be calculated as 

j ± = (^/2)F{±(a) ± (tt±/*)(C x a)), 



where 



(A 1 F 2 V/3TrN)(2m 3 )^ 2 h~ 3 E 1 k 



/2B± 



B± = ^- T f "{(ri ±r 2 J(3F 2 + 4F 2 cos6 
°™ Jo Jo 

+(r x ± r 2b )(-4:FoF 1 - 8F? cos (9)} sin 3 6d6d<f>. 



(4) 



and 



1 /*27r 
A = ^lo Jo {l J ( k ' k ')! 2 Pi( k '- k ) 

+ P(k, k')} (1 - cos 6) sin 6d6d(j). 

depends on all components of momentum transfer as 
it includes spin-spin correlations between different ions. 
and are the angles for spherical coordinates with re- 
spect to the longitudinal current, cr^ 1 is the resistivity 
given in terms of the integral A , 

cr" 1 = p\\ = (37rm/2he 2 )(V/N)(A/E F ). (5) 



n=0 
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Thus, the spin current j s 
j~ are 



j+ -j 



and the charge current 



a Inverse spin Hall effect near T 



j a = H Ta + <7,| + *_)/(2*)}(C x a), (6) 

f = a||J-{(* + -*_)/(2$)}(Cx<7). (7) 

As schematically shown in Fig. [3a), f° r the incident spin 
current in the ISHE configuration, the spin-up (+) and 
spin-down (— ) conduction electrons are scattered by skew 
scattering to the same side with the transition probability 
proportional to the terms i\ + r 2 (r 2 represents r 2a and 
r2b) and — n + T2, respectively. As a result, the Hall 
current j c is proportional to r 2 ■ 

The Hall resistivity is defined as pish = E±/jn. It has 
j\\ = a\\T from Eq. ©, and J5j. = - *_)/(2$) 

from Eq. ([7]). In the direction £ x cr, the spin current 
does not contribute to the electric field E±. Combining 
Eqs. dU and (0, the inverse spin Hall resistivity pish 
can be obtained as 



Pish 



KiB 2 F 2 (V/N) 2 (2m 5 ) 1 ' 2 /{2h:- 



e 2 4 /2 ), 



(8) 



where £?2 = {B+ — B-)/2. B 2 involves only r 2a and r 



2b- 



B 2 = — 



/ / {r 2a (3F* + 4F? cos ( 
Jo Jo 



8ttj Jo 
+r 2b (-4F Fi - 8F? cos (9)} sin 3 

We find that the temperature variation of the Hall re- 
sistivity pish for the inverse spin Hall effect in a ferro- 
magnetic metal is that of B 2 , i.e., fourth-order correla- 
tions. The other factors in Eq. (jHJl do not show a signifi- 
cant temperature dependence close to the magnetic phase 
transition [1] . Note the on-site term r 2a = r^ b is indepen- 
dent of angles 9 and <j) so the integral gives, in a purely 
local approximation, p° SH oc Kir 2 (F 2 — AF F 1 /3)F 2 . 
This is analogous to the local approximation p xy for the 
anomalous Hall resistance obtained by Kondo, but ap- 
pears to compare less well with experiment Q: For the 
ISHE the non-local terms are more significant. 

Spin Hall effect. For this effect [9j a transverse spin 
current is generated from a charge current, so we can 
deduce the current equations from those of the inverse 
effect by interchanging the spin and charge labels of the 
currents and replacing the spin diffusion field J- by the 
driving external electric field £. To define a spin Hall 
transport coefficient we must specify the geometry of the 
measurement. We consider the non-local device scheme 
where the potentials are measured for parallel (Vf) 
and antiparallel (V^ p ) magnetizations of the injecting 
and detecting ferromagnetic electrodes. Both Vf and 
Vf p are proportional to £(* + -'I , _)/(2$) with different 
coefficients of proportionality. \&_|_ and follow from 
Eqs. ([3]) and (JSJ) as before, so the spin accumulation sig- 
nal psh = {V± — V± P )/j\\ K Pish- The proportionality 
constant depends on the transparencies of the barriers, 
which should be essentially independent of temperature, 
so that the temperature variation of the spin Hall resis- 
tivity psh is just that of the inverse spin Hall coefficient 
Psh, i-e., B 2 . 
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b Anomalous Hall effect near T, 




FIG. 1. (Color online) Schematic picture of (a) ISHE and (b) 
AHE near Tc- The skew scattering amplitudes for spin- up 
and spin-down conduction electrons (represented by arrows on 
the deflecting arrows) are proportional to r\+r2 and — n+ra, 
respectively. The potentials measured depend on correlations 
of 4 (respectively 3), local spins in (a) [respectively (b)]. 



Anomalous Hall effect. We also generalized the expres- 
sion for the AHE to include short-range spin-spin correla- 
tions. By interchanging the spin and charge labels of the 
currents and replacing the spin diffusion field T by the 
external electric field £ = — with magnitude £ and 
direction a for the spin states ±, the current densities [|| 
can be calculated as 



t = (<r\\/2)£ <t+(*±/$){Cx<t} 



Thus, the spin and charge currents are obtained as 

j s =a||£{(* + -vI/_)/(2$)}(Cxa), (9) 
f = <j l{ £a + a ll £{(* + + *_)/(2<&)}(C x a). (10) 

As shown in Fig. QJb), for the incident charge current in 
the AHE configuration, the spin-up (+) and spin-down 
(— ) electrons are scattered by skew scattering to opposite 
sides with the transition probability proportional to the 
terms n +r 2 and — n +r 2 , respectively. As a result, the 
Hall current part in j c is proportional to r\. 

The Hall resistivity is defined as pn = E±/jn. Now 
j|| = a || J- from Eq. ©, and E± = - #_)/(2$) 

from Eq. (fTTTj) . In the direction £ x a, the spin current 
does not contribute to the electric field E±. With Eqs. 
(|4|) and ((5|), the anomalous Hall resistivity pn can be 
obtained as 

p H = A 1 (V/N) 2 (2m 5 ) 1 ^B 1 F 2 /(2h 4 e 2 E 1 F /2 ) 1 
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where B\ = (B + + B_)/2. B\ involves only v\: 

B 1 = —\ / rA'SFl - 4F F 1 - AF? cos 9} sin 3 9d9d(j). 
SttJo Jo 

If we consider only the site-diagonal n = term in 
r t = rl, B x = r?(T 2 - 4F F 1 /3), and p H will reduce 
to Kondo's result. We remark that the definitions of the 
local r° and r§ terms were given by Kondo. Only the 
r° term was found to contribute to the AHE Q. The 
effect of the r\ term, on the other hand, has been hid- 
den for about 50 years. Of course Kondo's paper was 
written almost a decade before Dyakonov and Perel's [l| , 
whose ideas have been explored even more recently with 
experimental manipulation of spin currents. 

The difference between the AHE and the ISHE near Tc 
originates from the different symmetries of the incident 
charge current in the AHE and incident spin current in 
the ISHE, as in Fig. 1. This symmetry difference gives 
rise to the distinct orders of spin-spin correlation: Near 
Tq the temperature variation of ph in the AHE is de- 
termined by a three-spin correlation n, while the Pish 
in the ISHE is determined by four-spin correlations r^a 
and r2b- The difference in order can be traced back to 
the different terms linear in spin orbit in Eq. ((T|): For 
charge current the term linear in local spin S„ combines 
with two exchange terms (for the second-order Born ap- 
proximation there are three powers of the Hamiltonian) 
to give a third-order term in local spins, whereas for spin 
current any of the last three terms, with two local spin 
operators, combine with two exchange terms to give a 
four-spin operator. In particular, the third-order corre- 
lations vanish by symmetry above Tc whereas the fourth 
order does not. In the limit of zero momentum transfer, 
i.e., k ~ k' = 0, ?'i and T2 a (=?"2&) simplify into first- (xi) 
and second-order (^2) non-linear uniform susceptibilities 
flTj . respectively, as shown in Fig. 2. The divergence at 
Xc should be smeared out in the r\ and r^a and rib due 
to the finite momentum transfer k — k', as for the resis- 
tance [l2j]. As a result, the different symmetries should 
be reflected by a single peak below Tc in the ph and two 
peaks of opposite sign above and below Tc in the Pish- 

Discussion. This theory is formulated for a pure crys- 
talline metal and it would normally be considered as an 
"intrinsic" mechanism. In fact, the dependence of the 
Hall constant on longitudinal resistance resembles more 
what is classified as "extrinsic" models for skew scat- 
tering from impurities. This is not an accident, since 
it is the fluctuations of the local moments of d orbitals 
which skew scatter close to the Curie temperature. Fert 
and co-workers (l3| . in a study of dilute alloys which 
are truly "extrinsic" in the sense of including impuri- 
ties, approached the transport in a way that looks dif- 
ferent, but is actually compatible. Rather than finding 
scattering amplitudes from microscopic potentials, they 
replaced the d orbitals by a set of phase shifts to be 
determined empirically. They included p- and d-wave 
phase shifts, more significant in the alloys, rather than 




FIG. 2. (Color online) Schematic behavior of the anomalous 
(left) and spin Hall resistances (right) around Tc- In each case 
the power-law divergence (dashed lines) of the corresponding 
non-linear susceptibility xi ° r X2 i s cu t on " by the finite Fermi 
surface. The different shapes (solid lines) of the anomalies 
reflect that the third-order correlations vanish about Tc and 
the fourth-order correlations change sign across Tc- 



the s- and p-scattering channels included here. As it is 
enough for skew scattering to include interference from 
phase shifts of different parities, the real difference is in 
the inclusion of resonant effects. We remark also that 
expression pn oc t°(Fq — 4FqFi/3)AiF2, obtained in the 
site diagonal limit is the equivalent of that of Fert et 
at ph oc Xd sin S p sin (25d — S p ) sin 2 Sd with A1F2 (com- 
ing from the scattering amplitude of the p-wave exchange 
amplitude and the local spin via the spin-orbit coupling) 
corresponding to the change in phase shifts for different 
components of d-wave scattering from the resonant po- 
tential \d sin 2 Sd- The factor (TJ 2 — 4ToTi/3) corresponds 
to the term sin(5 p sin (25d — S p ) 7 while the former term 
representing interference between the s wave (Fo) and p 
wave (Fi), and the latter term representing interference 
between the p wave and d wave. The similarity ends here: 
There is no equivalent of r° in Fort's theory, which was 
formulated for alloys at low temperatures, while here we 
are, as in Ref. Q, interested in the collective spin fluctua- 
tions near criticality. The distinction between the anoma- 
lous and spin Hall resistivities depends on the difference 
between the r 1 and the two r 2 terms. We remark that 
many-body effects can be included by use of the calcu- 
lated temperature-dependent phase shifts for an isolated 
quantum impurity [|| or in lattice models [l4[ for mixed- 
valence systems where there is a temperature-dependent 
effective phase shift. 

Our theory includes skew scattering but neither side- 
jump nor the "topological" lattice effects coming from 
anomalous velocities generated by spin-orbit terms ffol 
[Tr| and should thus be applicable when the conductivity 
is large [l7|. The band theories can be extended to fi- 
nite temperatures by replacing the magnetization of the 
band by the temperature-dependent value pn oc p 2 (Af) 
but it is hard to see how the anomalous velocity could 
produce non-monotonic behavior above Tc, in contrast 
to our approach. 

Comparison with experiments. Recently Wei et al. Q 
have performed an inverse spin Hall effect experiment in 
a weak ferromagnetic alloy near Tc , where a dip and 
a peak are observed in the spin Hall resistivity pgn be- 
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low and above the Tq, respectively. We argue that our 
theory can qualitatively explain the experimentally ob- 
served anomalous behaviors in psh near Tq. The theory 
we provide here including correlations between spins on 
more than one atom seems to be necessary to explain 
the experiments. From simulations Q on short-range 
Heisenberg models the on-site terms are rather smooth 
and only the off-site terms are non-monotonic. Thus fea- 
tures close to Tq may be dominated by the correlations 
that are off-diagonal in the site index. 

Conclusions. We have presented a theory of the in- 
verse spin and spin Hall effects that provides an essen- 
tial complement to the classic theory of Kondo for the 
anomalous Hall effect, which has remained unchallenged 
for 50 years. We have found an essential difference be- 
tween the anomalous Hall effect and the inverse spin Hall 
effect near Tq in a ferromagnetic metal. Our theory can 
be compared to recent experimental results for the SHE 



in ferromagnetic alloys. 
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